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ABSTRACT 


A model  is  developed  for  sliding  and  rolling  friction  that  can  be  used 
in  simulations  of  dynamic  systems  involving  mechanical  elements  that  are 
subject  to  friction.  Sound  bases  for  the  model  are  developed  around  the 
hypothesis  that  the  origin  of  friction  is  in  quasi  static  contact  bonds  that  are 
continuously  formed  and  subsequently  broken.  Sliding  friction,  including 
coulomb  friction  and  the  so-called  "stiction,"  as  well  as  rolling  friction  are 
believed  to  be  realistically  simulated  by  the  model.  The  model  has  been 
proven  to  accurately  simulate  transient  torque  and  rate  measurements  on 
ball  bearings  tested  in  the  laboratory. 
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I.  INTRODUCTION 


Recently,  a series  of  ball  bearing  tests  were  conducted  by  H.  Shibata 
in  the  Aerospace  Corporation's  Guidance  and  Control  Laboratory  under  the 
direction  of  D.  J.  Griep.  In  discussions,  Mr.  Griep  described  to  the  author 
in  detail  the  dynamical  behavior  of  a simple  experiment  with  three  balls  on  a 
flat  plate  with  another  plate  on  top  the  three  balls.  His  description  clearly 
indicated  that  small  displacements  of  the  balls  from  a rest  position  caused  an 
elastic  restoring  force  on  the  balls.  This  deduction  was  obvious  from  the 
observation  of  lightly  damped  lateral  oscillations  about  the  rest  position  when 
the  balls  were  disturbed.  Discussions  also  encompassed  the  stochastic 
nature  of  rolling  friction  which  the  Aerospace  tests  were  largely  concerned 
with. 

Mr.  Griep  welcomed  the  offer  of  this  author  to  attempt  to  model  and 
simulate  the  quasi  static  and  stochastic  friction  phenomena.  The  following 
report  is  the  result  of  this  relatively  meager  exposure  to  the  world  of  experi- 
mental friction.  In  the  search  for  a reasonably  plausible  simulation  model 
that  would  fit  the  experiment  observations,  some  latitude  has  been  taken  with 
fundamental  concepts  of  materials  properties  that  no  doubt  are  quantum 
mechanical  and  electrodynamical  in  nature  arid  so  the  readers  forbearance 
in  this  respect  is  solicited. 
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U.  A THEORY  OF  SOLID  FRICTION 

A solid  friction  simulation  model  is  desired  that  will  accurately  repre- 
sent the  physical  phenomena  commonly  referred  to  as  static  friction,  coulomb 
friction,  and  rolling  friction.  The  combined  effects  of  static  and  coulomb 
friction  are  usually  depicted  as  shown  in  Figure  1. 


Friction 
Force  of 
Torque 


t 


Stiction  Force 

— •♦——Coulomb  Force 


Velocity 


Figure  1 . Conventional  Representation  of 
Solid  Friction 

Although  generally  accepted  as  a basic  physical  macro-phenomenon, 
coulomb  friction  may  la  explained  on  the  basis  of  the  quasi  static  properties 
of  materials.  Typically,  these  strength  properties  of  solids  appear  ns  shown 
in  Figure  2. 

When  a solid  is  stressed  beyond  its  yield  point  the  material  undergoes 
permanent  deformation  and  when  the  stress  is  removed  the  material  behaves 
elastically;  i.e.,  stress  vs  strain  is  linear  during  stress  removal. 

Usually,  strain  is  thought  of  as  being  caused  by  an  imposed  stress. 
However,  the  reverse  appears  to  be  the  convenient  way  to  think  of  sliding 
friction  between  two  solid  bodies;  i.e.,  stress  or  force  or  torque  resisting 
the  initial  relative  motion  of  the  two  solid  bodies  is  caused  by  their  relative 
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Figure  2.  Strength  Properties  of  Solids 

motion  that  is  brought  about  by  an  overpowering  external  force  acting  on  the 
bodies.  Thus,  as  strain  or  relative  displacement  of  the  bodies  is  increased, 
an  internal  restoring  force  proportional  to  this  displacement  is  produced 
(considering  the  motion  between  the  two  bodies  to  be  such  that  "breakaway" 
has  not  occurred  so  that  the  two  bodies  behave  as  one  through  the  mechanism 
of  partial  cohesive  bonding).  Further  relative  displacement  causes  yielding 
or  inelastic  deformation  of  the  solid  materials  at  the  bonding  interface  while 
the  restoring  force  increases.  A maximum  restoring  force  is  reached  and, 
if  the  externally  applied  force  is  held  constant,  the  two  surfaces  start  to 
break  free  of  each  other  during  which  time  plastic  deformation  takes  place 
while  the  two  surfaces  are  still  in  solid  "tacky"  contact.  Finally,  rupture 
takes  place  and,  if  the  two  surfaces  are  not  physically  separated  by  some 
finite  distance,  the  contacting  elemental  surfaces  just  ruptured  will  bond 
and  break  continuously  with  other  elemental  surfaces  as  this  process  proceeds. 

This,  then,  is  a description  of  what  is  called  solid  friction.  It  pertains 
to  solids  in  contact  without  benefit  of  lubrication.  The  relationship  of  "stiction' 
and  coulomb  friction  to  this  solid  friction  now  becomes  clear.  The  coulomb 
force  is  related  to  the  interface  bond  rupture  stress  and  the  stiction  force  is 
related  to  the  maximum  or  iiltimate  stress  of  the  interface  bond.  Static 
friction  or  stiction  and  coulomb  friction  are  seen  to  be  indistinguishable  in 
materials  that  exhibit  brittle  fracture  properties  as  illustrated  in  Figure  3. 
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Figure  3.  Brittle  Material  Failure  Characteristic 

Ductile  materials  can  be  expected  to  exhibit  the  stiction  characteristic  shown 
in  Figure  4 where  the  "tacky"  phase  is  encountered,  which  causes  an  essen- 
tially unstable  progression  from  the  point  of  maximum  stress  to  the  point  of 
rupture.  It  is  evident  with  this  ductile  material  model,  then,  that  it  is 
possible  for  dynamic  relaxation  oscillations  to  exist  in  spring  mass  friction 
systems  due  to  the  negative  slope  portion  of  the  streSs-strain  curve.  It  is 
this  behavior  that  accounts  for  relaxation  oscillation  phenomena  such  as 
chatter,  squeak,  squeal,  and  stick-slip. 


Figure  4.  Ductile  Material  Failure  Characteristic 
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It  would  seem  too  much  to  expect  that  the  mechanism  described  above 
would  apply  equally  well  to  both  sliding  and  rolling  friction,  and  yet  it 
apparently  does.  Up  to  this  point  we  have  not  made  a distinction  between 
sliding  and  rolling  friction,  and  this  was  purposely  done  to  imbue  the  idea 
that  solid  friction  encompasses  both  of  these  types  of  friction.  The  difference 
between  these  types  of  friction  may  be  explained  simply  by  reference  to 
Figures  5 and  6. 

Figure  5 shows  the  case  of  sliding  friction  where  the  interface  cohesively 
bonded  regions  undergo  principally  shearing  stresses.  Thus,  the  interface 
bond  material  properties  that  govern  the  sliding  friction  process  are  mostly 
their  shear  properties. 


External  Force 


JT. 


7" 


Shear  stresses  in  interface  bonding  regions 


Figure  5.  Sliding  Friction 

On  the  other  hand,  the  case  of  rolling  friction  is  illustrated  in  Figure  6 
wherein  the  interface  cohesively  bonded  regions  are  compressively  stressed 
on  the  front  side  of  the  contact  area  and  are  stressed  in  tension  on  the  back 
side.  Thus,  for  rolling  friction,  tension  and  compression  properties  prin- 
cipally govern  the  friction  process. 

Sliding  friction,  then,  involves  shear  failure  because  it  is  the  only 
process  by  which  continuous  rupture  can  occur.  Tensile  failure  would 
separate  the  bodies  at  rupture  and  compressive  failure  would  break  one  or  both 
of  the  bodies.  Cutting  of  a solid  material  with  a tool  can  be  thought  of  as  this 
kind  of  process,  i.e.,  continuous  shearing  rupture. 
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Figure  6.  Rolling  Friction 

Rolling  friction  is  seen  to  also  involve  continuous  rupture,  but  the 
friction-associated  rupture  appears  to  be  tensile  in  most  engineering  appli- 
cations such  as  in  ball  bearings.  Rolling  mill  crushers  are  an  example  of 
the  continuous  rupture  process  that  involves  continuous  compressive  failure. 
Ball  bearings  no  doubt  fail  in  compression  when  overloaded,  but  during 
normal  operation  the  compressive  stresses  developed  on  the  front  side  of 
the  contact  area  merely  helps  establish  the  cohesive  bond  that  must  later  be 
broken  in  tension  on  the  back  side  of  the  contact  area. 

There  is  a question  whether  stiction,  per  se,  can  be  distinctly  identi- 
fiable in  rolling  friction.  The  author  believes  that  it  is  not.  This  preconcep- 
tion is  based  on  the  notion  that  the  interface  bond  strain  is  distributed  over 
the  area  of  contact  where  fracture  occurs  (i.e.  , the  back  side  where  tension 
failure  proceeds)  and  thus  an  averaging  of  the  strain  over  the  stress. strain 
curve  apparently  takes  place. 

The  phyrics  of  the  continuous  stress  rupture  process  is  not  understood 
but  conjecture  leads  this  writer  to  believe  that  an  explanation  perhaps  lies  in 
Dorn's  theory  (References  4 and  5)  that  creep,  leading  to  stress  rupture, 
occurs  in  thermally  activated  deformation  processes  or  perhaps  vice  versa, 
i.e.  , as  in  the  friction-caused  deformation  process  which  induces  thermal 
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activation  and  local  interface  heating.  Thus,  the  back  side  bond  creep  rates 
being  forced  by  rolling  velocity  are  the  result  of  contact  bond  energy  of 
activation  and  entropy  of  activation  produced  by  rolling  compression  on  the 
front  side.  The  entropy  of  activation  is  thus  believed  to  be  the  thermal 
energy  generated  by  the  rolling  friction  process. 
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m.  A SIMULATOR  FOR  SLIDING  FRICTION 
IN  DUCTILE  MATERIALS 

The  proceei  juet  described,  coupled  with  Newton's  second  law,  pro- 
vides a model  of  solid  sliding  friction  that  can  be  represented  by  the  sche- 
matic diagram  of  Figure  7.  This  process  is  simulated  by  the  block  diagram 
of  Figure  8 where  the  net  force  available  to  accelerate  the  mass  is  the  dif- 
ference between  the  component  of  the  external  force  parallel  to  the  contacting 
surface  elements  and  the  frictioiial  force.  The  time  integral  of  the  net  force 
produces  the  velocity  ic  which  is  integrated  to  find  the  elemental  displacement 
Ax  from  the  initial  equilibritun  position.  This  elemental  displacement  is 
redefined  each  time  the  direction  of  motion  changes.  Thus,  when  relay  R^ 
in  the  diagram  is  activated  by  positive  velocity,  contact  ia  closed,  and 

contact  R^  is  closed  prior  to  fracture,  which  permits  integration  to  obtain  the 
-tAx  relative  displacement  from  the  equilibrium  position.  When  Ax  reaches 


X 

F 

tang 


Figure  7.  Dynamical  Sliding  Friction  Model 
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Figure  8.  Sotid  Friction  Simulator 


fpj  (assuming  AXpj^  = 0),  relay  Rj  is  actuated,  corresponding  to  fracture  of 
the  contact  bond,  and  the  relay  contact  is  opened  thereby  holding  the  value 
of  Ax  on  the  integrator  output.  The  fracture  level  of  stress  is  then  filtered 
and  multiplied  by  the  contact  area  to  obtain  the  net  solid  frictional  force. 

This  force  is  obviously  constant  as  long  as  x is  positive  and  fracture  has 
occurred  and  is  seen  to  correspond  to  the  coulomb  friction  force. 

If  the  applied  force  is  reversed,  i will  change  sign  causing  closure  of 
the  R^^-  contact.  The  R^  contact  will  be  closed,  because  will  be  positive, 
corresponding  to  the  establishment  of  the  contact  bond  in  the  reverse  direc- 
tion. Thus,  Ax  will  start  to  decrease  from  the  value  Ax  = * ,• 

Pf 

For  X < 0,  the  ^■23'*'  contact  is  opened  and  the  ^33*  contact  closed  so 

that  the  friction  level  is  set  by  the  negative  stress  function  o-^.  At  the  instant 

of  X reversal,  r = a due  to  the  fact  that  or  is  slaved  to  o’  by  servo  action; 
n p n p ^ 

i.e. , is  compared  to  cr^  and  any  error  is  integrated  to  bias  the  input  to  the 

negative  stress  function  generator  so  as  to  drive  the  error  between  o’  and  o’ 

n p 

to  zero.  It  may  be  noted  that  the  positive  stress  function  o’^  is  slaved  to  o’^ 

for  negative  motions  (i  < 0)  when  o-^  is  being  read  out  to  generate 

Relay  switching  with  contacts  Sz  and  R^2  is  necessary  to  implement  the 

load  removal  characteristic  shown  in  Figure  2.  A filter  with  time  constant 

is  added  in  the  simulator  to  filter  the  R^^  switch  transients. 

The  net  friction  level  is  proportional  to  the  load  normal  to  the  surface 

in  the  case  of  static  and  coulomb  friction.  The  exact  mechanism  of  this 

relationship  is  not  understood  fully  and  so,  for  the  purposes  of  simulation, 

it  is  assumed  that  the  contact  or  bond  area  is  proportional  to  the  force  normal 

to  the  surface.  Thus,  in  the  block  diagram  represents  this  force. 

By  multiplying  the  stress  by  area  to  obtain  the  frictional  force,  it 

becomes  possible  to  introduce  friction  noise  in  what  is  believed  to  be  a 

realistic  manner.  The  area  noise,  A . , is  that  contact  area  variation 

noise 

produced  by  surface  irregularities.  It  is  a stochastic  function  of  the  variable 

X and  so  a random  function  generator  is  shown  that  generates  A (x)  for 

noise 

purpose. 
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Laminar  (viicoua)  and  turbulent  friction  terms  have  been  included 
in  the  simulator  block  diagram  for  generality.  Laminar  friction  is  assumed 
to  be  propor^onal  to  velocity,  and  turbulent  friction  is  assumed  to  be  pro- 
portional to  the  square  of  velocity. 


IV.  A SIMULATOR  FOR  SOLID  FRICTION 
IN  BRITTLE  MATERIALS 


An  alternate  and  much  simpler  fric^on  limulator  than  that  in  Figure  8 
is  shown  in  Figure  9.  This  simulator  should  be  applicable  to  rolling  friction 
and  brittle  dry  sliding  friction  and  perhaps  to  lubricated  sliding  friction 
involving  ductile  materials.  In  this  simulator,  advantage  is  taken  of  the 
observation  that  the  friction  stress  function,  f,  being  a function  of  displace- 
ment, X,  can  be  differentiated  with  respect  to  time 


df  df  dx 

3T  = ai  * ar  = 


f'i 


and  the  result  integrated  to  generate  f.  It  is  observed  that  V is  easily  found 
as  a function  of  x.  The  simulation  block  diagram,  however,  shows  f'  being 
generated  using  a itinction  generator  with  f as  the  input.  Thus  it  is  required 
to  determine  V as  a function  of  f.  This  determination  can  be  illustrated  with 
a few  examples . 

Example  1 

Assume  f(x)  = kx,  -x^  < x < x^ 

or  f(x)  = kx^,  )x|  > x^ 

Then  V - k,  -x^  < x < x^  , < f < kx^ 

= 0,  , |f  I > kx^ 
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Figure  9.  Friction  Simulator  For  Sliding  or  Rolling  Plus  Fluid  Friction 


This  is  easily  integrated  to  obtain 


X + c 


-TiTTTjr 


The  constant  of  integration  is  deterixiined  by  knowing  that  the  friction  f = 0 
ivhen  at  rest  at  x = 0.  Thus,  c = ^/Y  we  obtain 


(yf^jX  + 1) 


Referring  to  Figure  9,  the  time  derivative  of  the  friction  function  f is 
integrated  to  obtain  f(x).  Starting  at  x = 0,  f = 0,  and  for  it  > 0,  it  may  be 
seen  that  f(x)  asymptotically  approaches  f^  as  indicated  below.  When  x 
changes  sign  and  the  direction  of  motion  changes,  relay  is  actuated  thereby 
switching  R^^*  This  changes  the  sign  of  f from  -f  to  - causing  the  value  of 
f'  to  change  from  f'(f)  to  a value  f'(-f).  This  behavior  is  exhibited  for  all 
levels  of  friction. 


Reflecting  upon  the  example*  above,  one  see*  that  the  technique  of 
generating  f fail*  for  double -valued  function*,  provided  f'  i*  required  to  be 
continuou*;  i.e.,  the  integrator  become*  *tuck  a*  f'  pa**e*  through  zero. 
Thi*  might  be  circumvented  by  intentionally  placing  a diecontinuity  in  f ' 
where  it  pa**e*  through  zero  a*  i*  illuetrated  below: 


However,  the  eimulator  appear*  to  be  incapable  of  producing  double- 
valued function*  with  diode  function  generator*,  and  *o  only  monotonic  func- 
tion* can  be  eimulated  with  thi*  type  of  equipment.  It  is  believed  at  thi*  time 
that  monotonic  function*  will  adequately  represent  rolling  friction  but  not  the 
so  called  static  friction  associated  with  sliding. 

Since  monotonic  function*  appear  to  be  suited  to  simulation  by  the 
diagram  of  Figure  9,  the  characteristic  limiting  of  f when  f'  s 0 is  used  to 
great  advantage.  For  example,  in  rolling  friction,  it  can  continue  to  persist 
and  X will  increase  indefinitely  while  f is  limited. 


f: 


V.  SIMULATION  RESULTS 


A computer  aiikiulation  of  the  simple  model  shown  in  Figure  9 was 
implemented  on  an  EAI  TR-10  Analog  Computer.  The  computer  diagram 
used  is  shown  in  Figure  10.  Experimental  data  on  a 60  mm  bore  Class  7 
super  precision  12-degree  contact  angle  ball  bearing  with  a 60  lb  preload 
was  obtained  by  Hiro  Shibata  in  the  Aerospace  Guidance  and  Control  Laboratory 
under  the  direction  of  Dave  Griep.  The  experimental  apparatus  was  set 
up  to  drive  the  inner  race  Of  the  ball  bearing  with  a torque r whose  current 
was  closed-loop  servo-controlled  to  maintain  highly  accurate  angular  rates 
which  were  measured  'with  a gyro.  The  rate  could  be  commanded,  and 
several  tests  were  run  with  sinusoidal  input  rates.  Rearing  friction  torque 
was  measxired  with  a high  accuracy  Micro  Gee  position  table  adapted  for 
use  as  a torque  measuring  instrument.  The  measured  friction  torque  and 
gyro  rate  were  recorded  for  four  sinusoidal  inputs  to  the  rate  drive  servo. 

The  frequencies  of  these  four  runs  were  0.5,  0.1,  0.04,  and  0.01  cps  and 
the  angxilar  rate  amplitude  was  0. 125  deg/sec  in  all  four  of  these  runs.  The 
recorder  trace  data  was  reduced,  and  friction  torque  T^  was  plotted  versus 
angular  displacement  6 which  was  computed  using  gyro  data.  The  slopes  of 
the  resulting  curves  were  measured  to  obtain  dT^/d6  information  that  could 
be  plotted  versus  T^.  The  resulting  test  data  is  in  Figure  11.  Also  shown 
in  the  figure  is  the  function  approximated  by  the  DFG  in  the  TR-10  simulation. 

It  is  interesting  to  note  that  experimental  data  indicates  a square  law 
T'  vs  T relation.  An  empirical  curve  fit  equation  is  also  shown  in  Figure  11 
for  this  data.  The  empirical  equation  determined  is: 


wher'j 


dT 

3^ 


Y(T  - Tq)^ 


■YTq^  =■  36  in-oz/deg 
Tq  = 6 in-oz 
y =1.0  (in-oz-deg) 
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Figure  10.  EAI  TR-10  Analog  Computer  Diagram 


11.  Experimental  Data  and  Function  Fit  of 


Simulation  runs  were  made  employing  constant  rate  6 inputs  to  the 
simulator  and  other  input  functions  such  as  sawtooth  and  sine  functions. 
Because  it  had  been  expected  that  the  experimental  test  results  could  be 
interpreted  according  to  the  coulomb  friction  model,  the  data  obtained  from 
the  tests  were  plotted  as  friction  torque  versus  angxilar  rate.  Accordingly, 
several  runs  were  made  simulating  the  tests  and  the  computer  torque  and 
rate  variables  were  recorded  on  an  x-y  plotter. 

Restilts  of  the  simtilation  are  shown  in  Figure  12.  For  comparison, 
the  experimental  data  are  shown  for  half  of  a cycle  for  each  frequency. 

The  validity  of  the  model  presented  here  can  hardly  be  disputed  on  the 
basis  of  this  and  other  comparisons  that  have  been  made.  The  comparison 
of  Figure  12  shows  that  the  simulation  faithfully  reproduces  the  unexpected 
hysteresis  characteristic  first  noted  by  H.  Shibata.  It  is  apparent  by  exam- 
ination of  Figure  12  that,  as  frequency  is  lowered,  the  step-like  nature  of 
coulomb  friction  is  approached  with  the  width  of  the  hysteresis  loop  reduced. 
In  the  limit,  as  frequency  goes  to  zero,  the  step-like  discontinuity  is 
approached  and  the  hysteresis  loop  width  will  vanish.  Thus  we  see  that 
coulomb  friction  is  a special  case  of  the  more  general  solid  friction  model 
presented  here. 

It  is  interesting  to  speculate  that  the  friction  model  developed  here  is 
applicable  to  magnetic  and  other  forms  of  hysteresis.  Physical  processes 
that  exhibit  hysteresis  are  now  believed  to  be  characterized  by  the  functional 
dependence  of  the  ordinate  variable  on  the  time  integral  of  the  abscissa 
variable.  The  resulting  functions  behave  as  a brush  whose  bristles  m.;st 
bend  as  the  brush  moves  in  one  direction  and  then  flop  or  bend  in  the  opposite 
direction  if  the  motion  is  reversed. 


VI.  CONCLUSIONS 


A general  model  for  solid  friction  has  been  presented  that  has  been 
validated  by  ball  bearing  friction  experiments  data. 

The  heretofore  unexpected  hysteresis  characteristic  in  rolling  friction 
found  by  Shibata  and  Griep  was  somewhat  amazingly  simulated  by  the  simple 
computer  simulation  over  a wide  frequency  range. 

Coulomb  friction  was  found  to  be  a special  case  of  the  general  solid 
friction  model. 

Static  friction  (stiction)  is  exhibited  during  sliding  between  ductile 
materials. 

Stiction  does  not  occur  in  sliding  between  brittle  (or  well  lubricated) 
materials.  It  is  speculated  that  stiction  is  rare  or  completely  absent  in 
rolling  friction. 

Rolling  friction  involves  a compressing /tensing  process  and  sliding 
involves  a shearing  process. 

Solid  friction  is  caused  principally  by  contact  bond  stress  and  is  quasi 
static  in  character.  Making  and  breakin|  of  the  contact  bonds  is  dynamic  and 
governed  by  motion  of  the  bodies  driven  by  the  external  or  inertial  forces. 

The  friction  model  presented  is  believed  to  be  related  theoretically  to 
magnetic  hysteresis  and  the  simulator  described  can  obviously  be  used  to 
simulate  other  forms  of  hysteresis. 

The  stochastic  nature  of  various  forms  of  friction  is  thought  to  enter 
through  the  randomness  of  contact  area  variation  and  a seemingly  plausible 
way  of  including  this  in  the  simxilation  model  is  shown.  However,  the 
accuracy  of  this  representation  has  yet  to  be  validated.  Analog  computer 
random  function  generators  are  not  available  that  directly  produce  functions 
of  variables  other  than  time,  and  so  a digital  simulator  might  be  more  appli- 
cable for  this  purpose. 
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